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Abstract 

The analytic properties of the energy density p{t) of the cosmic 
fluid, and the Hubble parameter H{t), are investigated near to the 
future singularity t = tg assuming different forms for the equation of 
state. First, it is shown that the inclusion of quantum eff'ects coming 
from the conformal anomaly modifies the singularity. Thereafter, we 
consider the effect coming from a bulk viscosity in the fluid. The 
viscosity tends to reduce the magnitude of ts, but does not alter the 
singularity itself (the exponent). Main emphasis is laid on the simple 
case when the equation of state is p = wp, with w a constant. 



1 Introduction 

It is actually not more than about five years since the possible occurrence 
of a singularity of the universe in the far future - usually called the Big Rip 
- was first discussed in the literature [HOE]. To begin with, two different 
types of singularities were investigated. Later, it became clear that there are 
two possibilities more, so that present we know about four different types. 
For convenience, we quote from Ref. ^ the following classification: 

(i) Type I ("Big Rip"): For t ^ ts, a — > oo, p ^ oo, and \p\ — > oo, or p 
and p are finite at t = tg. 

(ii) Type II ("sudden"): For t — > t^, a ^ as, p ^ ps, and \p\ — > oo, 

(iii) Type III: For t ^ ts, a ^ as, p ^ oo, and \p\ — > oo, 

(iv) Type IV: For t ts, a —>■ as, p 0, \p\ — *• 0, or p and p are finite. 
Higher order derivatives of H diverge. 



^E-mail: iver.h.brevik(8)ntnu.no 

^On leave from Tomsk State Pedagogical University, Tomsk, Russia. Email: gorbimo- 
vaog@yandex.ru. 
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Here the notation is standard, a meaning the scale factor and tg referring 
to the instant of the singularity. These theories assume the cosmic fluid to 
be ideal, i.e. nonviscous. 

As often happens in physics, when the mathematical formalism encoun- 
ters singularities, this indicates that the description has been too simplistic. 
It becomes natural to ask: can the singularity be modified or softened by 
drawing other physical effects into consideration? The answer turns out to 
be affirmative, due to at least two reasons. The first is the effect of quantum 
corrections due to the conformal anomaly. It turn up in various cosmological 
applications: dynamical Casimir effect with conformal anomaly [8], or dark 
fluid with conformal anomaly [9]. Explicit calculation shows this anomaly 
to cause a softening of the singularity [SI |TU] . A second reason - the one to 
be treated in the present paper - is to take into account the bulk viscosity 
of the cosmic fluid. From a hydromechanical standpoint a generalization of 
the cosmic theory so as to encompass viscosity is most natural. As we will 
see, the viscosity also acts in the direction of softening. Viscous cosmology, 
beginning with the seminal paper of Misner [TT] seems generally to have at- 
tracted increased interest in recent years. Some references in this direction 
are [HI [la [la [ig [161 [m [18] . 

For convenience, and for reference purposes, we give in the next section 
a brief overview of the singularity application of conformal anomaly theory. 
Thereafter, from Section 3 onwards we embark on the core of the present 
paper, namely to include a bulk viscosity in the Einstein equations. 

2 Quantum effects, when p = —p — Ap^. 

We shall content ourselves by considering specific examples. Close to the 
borderline between the quintessence and phantom regions, the equation-of- 
state (EOS) parameter w = p/p is near to -1. One natural choice for the 
equation of state is the following [H [5l [6l [^ : 

P = -P-V, (1) 

where A and (3 are constants. Thus p always diverges when p becomes infinite. 
Note that this may be considered as a viscous dark energy. 

We shall distinguish between two cases, viz. 1/2 < (3 < 1, and (3 > 1. 
They may be classified as follows, referring to the notation in the previous 
section: 
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(l)l/2</5<l: The singularity is of Type I. The dominant energy 
condition (DEC), requiring p > 0, p ± p > 0, is broken. The EOS parameter 
approaches -1 for either sign of A, when p — > cxd. This case imphes that p 
diverges in a finite future. One obtains |5j 



p^^t,-t) — , H^{ts-t)^, (2) 

with ts the singularity time and H = a/a the Hubble parameter. 

(2) /5 > 1 : The singularity is of Type III. The DEC is also now broken. 
The magnitude of the EOS parameter goes to infinity, i.e., w +oo (— cxd) 
for A < (A > 0). Equations ^ still hold. 

We shall be concerned mainly with the case (3 > 1 m this section. 

Consider next the quantum contribution to the conformal anomaly. The 
complete energy density is ptot = p + Pa- Taking the trace of the conformal 
anomaly energy-momentum tensor, = —pA + "iPAi plus observing the 
energy conservation law, 

Pa + ^H{pa+Pa)=^. (3) 

we find that 

PA = -4:pA - Pa/H. (4) 
Thus we obtain for the conformal anomaly energy density [5] 

Pa = j dta'^HTA 

= _ J_ y dta^H[-12bH^ + 246'(-if2 + h^h + H^) 

- (46 + 66") {H + 7HH + + 12H^H)] , (5) 

where b, b' and b" are constants, occurring in the expression for the conformal 
trace anomaly, 

Ta = b{F + -nR) + b'G + b"R. (6) 

3 

Here F is the squared Weyl tensor and G the Gauss-Bonnet invariant. Ex- 
plicitly, if there are N scalars, N112 spinors, Ni vectors, N2 gravitons, and 
Nhd higher derivative conformal scalars, one has for b and b' the expressions 

AT + QNi,2 + 12N^ + QIIN2 - 8Nhd 
120(47r)2 ' 



3 



_ N + llN,/2 + Q2Ni + UIIN2 - 28Nhd 
360(47r)2 ' 

whereas b" is an arbitrary constant whose value depends on the regulariza- 
tion. 

The quantum corrected FRW equation is, with = SuG, 

^,H' = p + pA. (8) 

As R is large near t = ts, we assume {?>/ k^)H'^ <^ \pa\- Hence, p ~ —pA- 
In the presence of quantum effects, suppose that 

p = Po{ts-ty (9) 

as in Ref. [5j. From the energy conservation equation 

~ 1-/3 

ff - _P (f _ 1^1+7(1-/3) 

we get, by taking the most singular term, p ~ —(46 + 6b")HH. Integrating 
with respect to t we get p to be of the form Qj, with 7 = 4/(1 — 2/3). Thus 

4 3-2/3 
pr^{t,-t)^, Hr^{ts-t)^ (10) 

near t = tg. If we compare this with the case of no quantum corrections, we 
see that the energy density diverges more rapidly in the present case, whereas 
H becomes less singular. That means, the quantum corrections moderate the 
singularity. Note that 

(t^_t)T^^O, (11) 

P 

showing that our approximation neglecting the left hand side of Eq. ([8]) is 
consistent. 

The case /3 > 1 as assumed here, corresponds to singularity of Type III. 
The same properties may be found for the Big Rip (Type I) singularity. 



3 Viscosity-generated singularity, when p = 

-P - f{p) 

Assume now that the cosmic fluid is viscous, with a bulk viscosity (. In gen- 
eral, C = C(p)- Hereafter we leave out the conformal anomaly, so that p^ = 0. 
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The FRW equation becomes accordingly (3/K^)iJ^ = p, where the matter en- 
ergy density p incorporates the modification due to the viscosity. We assume 
the cosmic fiuid to be thermodynamically a simple, one-component, fiuid. 
We start by assuming the general form p = —p — f{p), with /(p) a general 
function. 

This case was analyzed in Ref. [IS]. The essential difference from the 
foregoing is that the energy conservation equation is now 

p + 3H{p + p) = 9CH^ (12) 

(note that the scalar expansion is ^ = 3H). We let t = denote present time, 
and let subscript refer to present time quantities. The viscosity-generated 
singularity time is denoted as tgc;. From Eq. (9) in [19] we have 

(13) 



Singularities of Type I or III occur if t = t^c corresponds to p ^ oo. 

We assume henceforth that C(p) C is a constant, and consider two 
specific choices for the function /(p). 

4 The case when f{p) = ap 

We consider this case first because it is simplest. It corresponds to an equa- 
tion of state in the form 

p = wp = -{1 + a)p, (14) 
with a a constant. From Eq. (17) in [19] we have 

ZJ pt/tc 



where is a 'viscosity time' 



tc = il'^'C)-'. (16) 



It is convenient to introduce two new nondimensional quantities x and K, 

^ = — , i^=^— • (17) 
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From Eq. fllSp it is seen that H becomes infinite when 



X = \n{l + Kx). (18) 

This transcendental equation for x can have two solutions. 

(i) Nonviscous solution. There is always the solution x = 0, corresponding 
to tc = oo or C = 0. This is accordingly the solution for a nonviscous fluid. 
We will call the singularity time in this case tso- We get 

SaHo 

It is desirable to write the solutions such that the difference from the singu- 
larity time occurs explicitly. First, we have 

H(t) = ^ (20) 



which corresponds to 



,2/3a 



Next, as H oc p according to the FRW equation, we get 



{tso - ty 



P(i) = 77^^- (22) 



Comparing Eqs. fl20|) and fl22l) with Eq. ([2]), we see that corresponding equa- 
tions for H{t) and p{t) are of the same form, if /5 = 1. This is as we might ex- 
pect, as the equation of state in the present case corresponds to A = a, P = 1. 
The singularity is of Type I if a > 0, and of type III if a < 0. 

(ii) Viscous solution. In the viscous case, x > 1, Eq. (ITSl) has a second 
root if the slope K is larger than one. As K is actually the same as the ratio 
tso/tsc we thus get the following condition for solution 

is( < iso- (23) 

The presence of a bulk viscosity thus reduces the future singularity time. It 
would not be so easy to foresee this property beforehand. 

For a given value of tc, the transcendental equation ffTSl) does not allow 
us to solve for t^c explicitly. However, we may easily handle analytically the 
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limiting case t — t^,^, which is the situation of main interest here. Using the 
approximation e"!^" = 1 + r/tc with t = tsc^ — t, we obtain 

H{t)^^ (24) 

This formula holds for all values of (. 

It is thus apparent that we get the same singular behavior of H{t) in the 
viscous case as we did previously, in Eqs. fl20l) and ([2]) (the latter under the 
assumption that /5 = 1). Correspondingly, we get in the present case 

a(t) ~ (t,^-t)-2/3a^ ^^^^^^ (25) 

p{t) ^ {t,^ - t)-\ t^t,^. (26) 

Whereas the viscosity tends to shorten the future singularity time, it does 
not modify the exponents in the singularity. Again, the singularity is of Type 
I if a > 0, and of type III if a < 0. 

When dealing with viscosity we have introduced two different times, 
namely a viscosity time tc defined in Eq. f[T^ . and a singularity time t^c- 
One may wonder: can anything be said about the relative magnitude of 
these times? As most cosmological theories are formulated without viscosity, 
it is natural to assume that the cases of main physical interest pertain to 
setting ( very small, or tc very large. Accordingly, the quantity x defined in 
Eq. (IT7|) should be expected to be small, corresponding to -C tc- 

Another related point is to inquire whether there is a relationship between 
the singularity time tg with or without conformal anomaly considered in 
Section 2, and the viscosity singularity time tg^. Although there does not 
seem to be a close relationship, it is in our opinion most natural to identify 
the nonviscous time tgo in Eq. flTIJl) with the previous time ts- As we have 
seen, the effect from bulk viscosity is to reduce the singularity time. We thus 
suggest that the relationship 

tsC < ts (27) 

is quite general. 

Note that one may take into account quantum effects in this case too. 
However, the corresponding equations are very cumbersome and require nu- 
merical study. Such a study indicates that that quantum effects may again 
soften the future singularity. 
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5 Remarks on the case when f{p) = Ap^ 



In this more general case we go back to Eq. f|T3l) (still assuming that ( is 
constant). Leaving at first the integration limits unspecified, we have 



1 1 
3^C 



dp 



P 



1 + 



P' 



,/3-l/2 



(28) 



As (3 > 1/2, it is seen that the last term in the denominator (the viscosity 
term) dominates for large p. We perform the integration from an initial 
value of p already lying in this region, up to p ^ oo. The corresponding 
time integration limits are t and t^^. We obtain approximatively 



ts( — t 



1 1 



V3kA(3-1/2 p/3-1/2- 



Thus 



pit)r^itsC-t)- 



which generalizes Eq. fl26|) and reduces to it in the case when (3 
( !30l) has actually the same form as Eq. ([2]). 



t 



(29) 
(30) 

1. Equation 



6 Summary 

We may summarize our results as follows. 



1) We started out by giving a brief overview of the four types of future 
singularities, listing in Eq. ([2]) the known forms for p{t) and H{t) in the case 
where the equation of state is assumed as in Eq. ([T]). Actual in our context 
were singularities of Type I and Type III. 

2) Inclusion of quantum effects coming from the conformal anomaly - cf. 
Eq. ([5]) - caused the exponents of the future singularity to be modified. 

3) In Section 3 we left out quantum effects, but considered instead the 
influence from a bulk viscosity Assuming the simple form ( JT^ for the 
equation of state, we showed in Eqs. fl2^ - fl26|) that the viscosity acts so 
as to shorten the singularity time, but it does not change the nature of 
the singularity. The same property was found when dealing with a more 
complicated equation of state in Section 5. The explicit calculations were 
limited to the case when ( was a constant. 



8 



Finally we mention that, under quite general conditions, preliminary nu- 
merical work on quantum effects indicates that they tend to soften the future 
singularity. 
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